r ~i 

PHYSICAL PROPERTIES OF 
STRANGELETS 



Jes Madsen 

Institute of Physics and Astronomy, University of Aarhus, DK-8000 Arhus C, 

I >• II Ill-Ill. . /. «;» 1 1 < 1 1 1 . < I < 1 1 1 . < I It 



Low-baryon number lumps of strange quark matter, strangelets, are 
presently searched for in ultra-relativistic heavy-ion experiments at 
CERN and Brookhaven. This paper gives an overview of the physi- 
cal properties of strangelets with emphasis on experimental signatures 
such as mass and charge. Direct solutions of the Dirac equation with 
MIT bag model boundary conditions are applied and compared to cal- 
culations based on a smoothed density of states, as well as to simple 
mass-formulae based on the MIT bag. Most strangelet properties can 
be understood within an approximation including bulk, surface and 
curvature energy contributions. 

in 

X> INTRODUCTION 

A range of strong interaction parameters exists for which bulk strange quark 
£" — matter (SQM) (1,2) is absolutely stable relative to nuclear matter (E/A < 
930MeV), or metastable on a weak interaction time-scale. The reason is 
that a third Fermi-sea, that of strange quarks, makes it possible to lower the 
energy relative to systems composed of up and down quarks only. To create a 
more stable system, the energy gained must first compensate for the current 
mass of the strange quark, but since the typical Fermi-energies involved are 
w nuc i eon /3 ?s 310MeV, this can be the case if the strange quark is not too 
_/ massive (the current mass of the strange quark is believed to be 100-300 
^MeV). 

f^h (Meta)stability of bulk SQM has far-reaching consequences for neutron star 
*p^physics, and perhaps cosmology (for a recent review with emphasis on the as- 
O ^rophvsics applications, see (3); for a general overview on SQM, including an 
1} extensive list of references up to mid-1991, see (4); for updates on some SQM- 
related issues, see (5); an overview of strangelet experiments is given in (6)). 
Attempts to produce SQM in the laboratory necessarily focus on low-baryon 
number objects (strangelets). These are in general significantly destabilized 
relative to bulk SQM due to finite size effects, so even if strangelets do not ap- 
pear in the ultra-relativistic heavy-ion experiments at CERN and Brookhaven, 
this does not rule out the (meta)stability of SQM in bulk. Fortunately most 
of the experiments are sensitive to strangelets even if they decay on a weak 
interaction time-scale, and the experimental signature (first of all a very low 
charge-to-mass ratio) is quite unambiguous (though confusion with strange 
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hadronic matter, discussed by Dover elsewhere in this volume, is a possibility), 
so the next few years will very likely either prove the existence of strangelets 
or put significant constraints on strong interaction model parameters. 

In the following I will introduce a physical framework for understanding the 
main properties of strangelets. Starting with a brief discussion of the MIT 
bag model and shell model calculations for strangelets, I will go on to describe 
a simple strangelet mass-formula, which I later derive explicitly from the MIT 
bag model. This allows a comparison of strangelet properties as a function 
of parameters. Finally I stress the consequences of the results in relation to 
experimental efforts. The main message is that some strangelet properties 
(like low charge-to-mass ratio) are quite robust, but that for instance the 
masses (and therefore also the stability properties) are very uncertain due 
to the freedom in choosing parameters. Some parts of the manuscript are 
updated versions of another recent review (7). 

STRANGE QUARK MATTER WITHIN THE MIT BAG MODEL 

In the simplest version of the MIT bag model (8) noninteracting quarks are 
confined in a spherical cavity of radius R. They satisfy the free Dirac equa- 
tion inside the cavity and obey a boundary condition at the surface, which 
corresponds to no quark current flow across the surface. The bag itself has an 
energy of BV , where the bag constant, B, is a measure of the false vacuum 
energy that confines the quarks (B thus in practice acts like an external pres- 
sure on the bag). In the simplest version the energy (mass) of the system is 
given by the sum of the bag energy and the energies of individual quarks, 



Here the momentum h K> i is found from the Dirac equation. For states with 
quantum numbers (j, /), k = ±(j + -|) for / = j ± For a given quark flavor 
each level has a degeneracy of N Kt i = 3(2j + l) (the factor 3 from color degrees 
of freedom). For example, the 151/2 ground-state (j = 1/2, / = 0, k = — 1) 
has a degeneracy of 6. 

The level filling scheme is somewhat cumbersome. To find the ground state 
strangelet for fixed bag constant, quark-masses, and baryon number, one must 
fill up the lowest energy levels for a choice of radius; then vary the radius until 
a minimum energy is found (dE/dR = 0). Since levels cross, the order of levels 
is changing as a function of R. 

Results of such calculations are shown in Figure 1. Notice that the energy 
per baryon smoothly approaches a bulk limit for A — ► oo, whereas the energy 
grows significantly for low A. For low s-quark mass a very significant shell 
appears for A = 6 (3 colors and 2 spin orientations per flavor), and less 
conspicuous ones for a range of higher baryon numbers. As m s increases it 
becomes more and more favorable to use u and d rather than s-quarks, and 
the "magic numbers" change; for instance the first closed shell is seen for 
A = 4 rather than 6. 




(1) 
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FIG. 1. Energy per baryon as a function of baryon number for B 1 ^ = 145MeV 
and s-quark mass from 50-300 MeV in steps of 50 MeV (higher curves have higher 
m s ). Spiky curves are based on direct mode-filling calculations, smooth curves on 
the smoothed density of states from the multiple reflection expansion. The smooth 
curves are thus not fits to the shell model calculations, but derived from a smoothing 
procedure within the MIT bag model. 
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Equation (1) can be modified by inclusion of Coulomb energy, zero-point 
fluctuation energy, and gluon exchange energy. The zero-point energy is 
normally included as a phenomenological term of the form —Zo/R, where 
fits to light hadron spectra (8) indicate the choice Zo = 1.84 for = 
5 1/,4 /145MeV = 1. In the MIT bag model the zero point energy is to some 
extent a fudge factor adjusted to fit the spectra. Roughly half of this phe- 
nomenological term has a physical basis in a center-of-mass motion, which 
can be subtracted in a more direct fashion. The proper choice of Zo (or even 
the functional form of the zero-point energy) is not straightforward. As dis- 
cussed by Farhi and Jaffe (9) the value is intimately coupled to B and m s , 
as well as to the strong fine-structure constant, a s , and it is not obvious that 
values deduced from bag model fits to ordinary hadrons are to be preferred. 
For reasonable parameter values one sees a significant effect of the zero-point 
energy for A < 10, but the term quickly becomes negligible for increasing A, 
where its contribution to the energy per baryon goes like A~ A I 3 . Because of 
the great uncertainty, no zero-point energy is included in the numerical results 
presented in this overview. 

Bulk systems are described in terms of 3 parameters: The bag constant, B; 
the strange quark mass, m s (up and down quarks are normally assumed to be 
massless); and a s , the strong fine-structure constant, which describes gluon 
exchange interactions. For most purposes a non-zero a s can be "absorbed" 
in a reduction of B, so in the following I shall concentrate on a s = 0. A lower 
limit on B (5 1 / 4 > 145MeV) can be inferred from the stability of ordinary 
nuclei relative to up-down quark matter (ordinary nuclei do not spontaneously 
decay into strangelets since this would require a high order weak interaction 
to make sufficient numbers of strange quarks). A bag constant smaller than 
(164MeV) 4 permits stable bulk SQM for sufficiently low m s , whereas the 
metastability window relative to a gas of A's goes to (195MeV) 4 . 



STRANGELET MASS-FORMULAE 



Ordinary nuclei are often discussed in terms of mass-formulae, where the 
energy per baryon is a sum of constituent masses, Coulomb energy, surface 
energy, symmetry energy, and perhaps other, less important terms. In a rather 
similar fashion, the energy per baryon for strangelets, E/A = e, can be written 
as a bulk term, e°, plus a number of finite size corrections, 

^ — ^ ~\~ ^Coulomb ~\~ ^surface ~\~ ^curvature ~\~ ^zero--- (2) 

In contrast to nuclei, fcouiomb is negligible for most strangelets due to the 
cancellation of u, d, and s quark charges. The Coulomb-term is important 
for determining the ground state composition of strangelets, but the energy 
itself is negligible in the mass-formula. The surface tension is important as in 
nuclear physics, but more important for most choices of bag model parameters 
is the curvature energy, which plays a negligible role in nuclei. Both terms 
are a consequence of the bag model boundary conditions, which result in a 
depletion of massive s-quarks near the surface (this is further discussed below). 
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In a first approximation strangelets have constant density (slightly higher 
than the nuclear matter density) for fixed bag parameters, so that baryon 
number A oc R 3 , E SUT f oc R 2 , and E CUTV cx R, where R is the strangelet 
radius. In this approximation (where we neglect e zeTO for the time being) 

e«e° + c surf A- 1 / 3 + c curv A- 2 / 3 , (3) 

with c sur f pa lOOMeV and c curv ps 300MeV as typical values. (The coefficients 
will be derived in the following subsections). 

Stability relative to a gas of neutrons requires e < m n , which may be written 
as A > A^ D , where 



abs 
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.abs ( Csurf + [c s 2 urf + 4c curv (m n - e )] 1 / 2 > 
mm V 2(m n - e°) J ■ ( } 

Stability at baryon number 30 requires a bulk binding energy in excess 
of 65 MeV, which is barely within reach for m s > lOOMeV if, at the same 
time, i/G?-quark matter shall be unstable. The proposed cosmic ray strangelet- 
candidates with baryon number 370 (10) would for stability require a bulk 
binding energy per baryon exceeding 20 MeV to overcome the combined cur- 
vature and surface energies. Absolute stability relative to a gas of 56 Fe cor- 
responds to furthermore using 930 MeV instead of m n , whereas stability rel- 
ative to a gas of A-particles (the ultimate limit for formation of short-lived 
strangelets) would correspond to substitution of m\ = 1116MeV. 

One can also calculate the minimum baryon number for which long-lived 
metastability with respect to neutron emission is possible. This requires 
dE CUTV /dA + dE SUT f/dA < m n — e°, or 

„meta _ ( C surf + [tjrrf + 3c curv (m n - e )] 1 ' 2 ^ 

mm " I 3(m n - e°) J ■ (5) 

To have A™f^ a < 30 requires m n — e° > 30MeV, which is possible, but only 
for a narrow range of parameters. 

However, it is important to notice, that shell effects can have a stabilizing 
influence. As stressed by Gilson and Jaffe (11) the fact that the slope of E j A 
versus A becomes very steep near magic numbers can lead to strangelets that 
are metastable (stable against single baryon emission) even for e° > 930MeV. 
Also, the time-scales for energetically allowed decays have not been calculated. 
Pauli-blocking is known to delay weak quark conversion in strangelets, and 
this will probably have a significant influence on the lifetimes. 



Smoothed density of states 

I will now show how the strangelet mass-formula can be explicitly derived 
within the MIT bag model. 

Strangelets with A <C 10 7 are smaller than the electron Compton wave- 
length, and electrons are therefore mainly localized outside the quark phase 
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(9), which means that they can safely be neglected in mass calculations. 
Strangelets therefore do not obey a requirement of local charge neutrality, 
as in the case of SQM in bulk. This leads to a small Coulomb energy, which is 
rather negligible for the mass-formula (less than a few MeV per baryon), but 
which is decisive for the quark composition and therefore the charge-to-mass 
ratio, Z/A, of the strangelet. A characteristic of strangelets, which is perhaps 
the best experimental signature, is that this ratio is very small compared to 
ordinary nuclei. Indeed, for m s = the ground state strangelet has equal 
numbers of all three quark flavors and is therefore charge neutral. Whereas 
Coulomb effects have been consistently included (12), I will leave out those 
terms in the equations below. 

In the ideal Fermi-gas approximation the energy of a system composed of 
quark flavors i is given by 

E = ^(fi,- + N iVi) + BV (6) 

i 

Here N{ and Hi denote thermodynamic potentials, total number of quarks, 
and chemical potentials, respectively. B is the bag constant, V is the bag 
volume. 

The proper calculation of the energy by an explicit solution of the Dirac 
equation with MIT boundary conditions was described earlier. A physically 
more illuminating approach, first applied to strangelets (without the curvature 
term, the significance of which was realized in (13)) by Farhi & Jaffe (9) and 
Berger & Jaffe (14), is to use a smoothed density of states rather than the 
explicit mode-filling scheme. In the multiple reflection expansion framework 
of Balian and Bloch (15), such a smoothed density of states can be written as 

tMS**^)-*^)'*-}' « 

where the area S = § dS = 4ttR 2 for a sphere, and curvature C = 

§ ^-^r- + -j^J dS = 87ri? for a sphere. Curvature radii are denoted R\ and 

i?2- For a spherical system Ri = R 2 = R. The functions fs and fc are given 
below. 

In terms of volume-, surface-, and curvature-densities, riiy, and n 8) c 

the number of quarks of flavor i is 

N i= -n 1 dk = n itV V + n it sS + n itC C, (8) 

Jo ak 

with Fermi momentum kpi = — mf ) 1 ! 2 = — A?) 1 / 2 ; A 8 - = mi/fr. 
The corresponding thermodynamic potentials are related by 

Qi=Qi iV V + Qi i sS + Q ii cC, (9) 

where d£li/dni = — iV» , and d^lij/d^i = —nij. The universal volume terms 
are given by 
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«.-,v = ((1 - A?) 1/2 (l - f A?) + §A? In 1 + (1 " A ' )1/2 ) , (10) 

n, y = 4(l-A 8 2 ) 3/2 . (11) 
The surface and curvature contributions are derived from (14) 



and (16) 




The explicit forms of 5 and rii s were first given in (14), and Sl 8) c and 
rii c in (16). The cumbersome formulae will not be quoted here, but their 
properties are shown in Figure 2. Note in particular that whereas 5 and 
&i,c ~^ for m —> n (no quarks of flavor i are present for m > n), and 
&i,s ~^ also for m —> 0, the curvature energy given by Sl 8) c is non-zero for 
massless quarks. This is very important because it means, that the nearly 
massless u and d quarks, while not contributing to the surface tension, give 
a very large contribution to the curvature energy. This effect was forgotten 
in strangelet studies prior to 1993, leading to mass-formulae that were too 
optimistic (significantly underestimating the energy per baryon) at the low 
baryon numbers accessible in ultra-relativistic heavy-ion colliders (13). 

An equilibrium strangelet is found by minimizing the energy, i.e. from 

= dE = (Vi,vdV + n iiS dS + n iiC dC + HidNi) + BdV. (14) 

i 

Minimizing at fixed N{ for a sphere gives the pressure equilibrium constraint 

B = - a,v - 1 £ ^ - E toi.c (15) 

i i i 

The "optimal" composition is found by minimizing the energy with respect 
to Ni at fixed A and radius, giving 

= J2^dNi. (16) 

i 

Together with Eq. (6) these constraints give the properties of a spherical quark 
lump in its ground state. 

The solution is compared with the shell model calculations in Figure 1. 
The agreement is remarkably good, indicating that all the important physics 
(apart from the shells) can be understood in terms of the surface and curvature 
contributions. 
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200 

m (MeV) 

FIG. 2. Surface tension and curvature coefficient as a function of quark mass for 

1 /3 

quark chemical potential of 300MeV. The upper curve is fi 8 s , while the lower curve 

is Ifi^d 1 ^ 2 multiplied by the sign of the curvature term. Numbers on both axes 
scale in proportion to fi. Note the significant curvature contribution for massless 
quarks. 
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Bulk approximation 

Simpler, approximate mass-formulae can be derived by using a bulk ap- 
proximation to the chemical potentials. This was the approach taken in the 
first detailed investigation of the strangelet mass-formula within the MIT bag 
model performed by Berger and Jaffe (14). They included Coulomb correc- 
tions and surface tension effects stemming from the depletion in the surface 
density of states due to the mass of the strange quark. Both effects were 
treated as perturbations added to a bulk solution with the surface contribu- 
tion derived from Eq. (12). 

While showing some of the important physics, this approach (apart from 
Coulomb corrections) predicted constant E/A versus A for m s — ► and m s — ► 
Us, in striking contrast to the shell model results, because the very important 
curvature term was not included. Below I discuss the bulk approximation 
with curvature. 

Including bulk terms only, the energy minimization, Eq. (15) (with A = A s , 
and superscript denoting bulk values), changes to 



B 



E 

i=u,d 

3 



1 (A*?) 4 



4tt 2 



A In- 



o\4 r 



4tt 2 
(1 - A 2 ) 1 / 2 



(i-a 2 ) 1/2 (i-|a 2 ) 



A 



(17) 



and the baryon number density is now given by 



E 



0\3 



+ i^(l_A 2 ) 3 / 2 



(18) 



A bulk radius per baryon is defined by 

R° = W^) 1 ' 3 . 



(19) 



In bulk equilibrium the chemical potentials of the three quark flavors are 



equal, fi° 



u° 



e°/3, where e° is the bulk energy per baryon. One 



may approximate the energy per baryon of small strangelets as a sum of bulk, 
surface and curvature terms, using the chemical potential calculated in bulk: 



nlsS + A' 



"TV. 



(20) 



where 5° = 47r(i?°) 2 J 4 2 / 3 and C° = 87r(i?°) J 4 1 / 3 . Examples for 5 1 / 4 
145MeV are (with s-quark mass given in parenthesis; all energies in MeV) 

e(0) = 829 + 0A~ 1/3 + 351A" 2 / 3 
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e(50) = 835 + 61A~ 1/3 + 277A~ 2/3 
e(100) = 852 + 83A" 1 / 3 + 241A" 2 / 3 
e(150) = 874 + 77A" 1 / 3 + 232A" 2 / 3 
e(200) = 896 + 53A -1 / 3 + 242A -2 / 3 
e(250) = 911 + 22A" 1 / 3 + 266A" 2 / 3 

e (300) = 917 + 0A- 1/3 + 295 J 4- 2/3 . (21) 

The bulk approximations above generally undershoot the correct solution 
with properly smoothed density of states (smooth curves in Figure 1) by 2MeV 
for A > 100, 5MeV for A « 50, lOMeV for A « 10 and 20MeV for A « 5. 
This is because the actual chemical potentials of the quarks increase when A 
decreases, whereas the bulk approximations use constant \i. For massless s- 
quarks the expression for e(0) scales simply as B 1 / 4 . The same scaling applies 
for m s > e°/3, where no s-quarks are present; in the example above the 
scaling can be applied to e(300). For intermediate s-quark masses both e and 
m s should be multiplied by B\l£ to scale the results. For B 1 / 4 = 165MeV one 
finds e(150) = 985 + 93A- 1 / 3 +265 J 4- 2 / 3 ; e(250) = 1027+46A- 1 / 3 +284 J 4- 2 / 3 . 
Coulomb effects were not included here. They have no influence for m s — ► 0, 
but change the results by a few MeV for large m s . 

In connection with the shell-model calculations I described the effects of a 
zero-point energy of the form —Zo/R, and claimed that it was important only 
for A < 10. This can be understood in the bulk approximation of constant 
H, because the zero-point term per baryon is proportional to A~ A I 3 compared 
to A -1 ! 3 and A~ 2 I 3 for surface and curvature energies. The full term to be 
added to the bulk approximation expressions for a given e° is: 

£zero = -Z (4/243tt) 1 / 3 [2 + [1 - (3m s /e )] 3 / 2 ] ^ e^" 4 / 3 , (22) 

typically of order -200ZoMeVA- 4 / 3 . 

EXCITED STATES 

The calculations discussed so far have focussed on finding the ground state 
properties of strangelets, i.e. the lowest mass state for a given baryon number, 
A. Figure 3, which is based on direct mode-filling including Coulomb energy, 
instead fixes A = 20 and displays contours of equal E/A for the whole range of 
possible strangelets with that baryon number (a total of 1891 states). What 
is quite interesting from an experimental point of view is the close spacing 
of states around that of minimum energy. More than a hundred different 
states are within lOMeV per baryon from the ground state. The ground state 
itself has slightly positive charge for all choices of parameters tried so far, but 
neutral as well as negatively charged states are very close in energy. Many of 
these states will be stable against strong decays, and since some weak decays 
are suppressed by Pauli-blocking, presumably many different configurations 
could be sufficiently long-lived to reach the detectors. This may be good 
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FIG. 3. Contours of energy per baryon in MeV for all 1891 strangelet configurations 
with A = 20, m s = lOOMeV, and B\H = 1. Contours are plotted as a function of 
charge (Z = N u — A) and the number of strange quarks, N s . 



from a production point of view, but also means that there is not going to be 
a single well-defined strangelet signature, but rather numerous (meta)stable 
states and many more possible decay modes! 

Related results are shown in Figure 4 where the mass-formula in the bulk 
approximation has been used to compare E/A for strangelets of equal mass 
and charge. As is seen there, even for fixed B, m s , A, and Z, there is a huge 
span in strangelet mass! High resolution mass measurements may therefore 
be of significant interest for comparison between theory and experiment. 

CONCLUSIONS 



The physical properties of strangelets can be summarized as follows: 

1) (Meta)stable strangelets are possible for a wide range of parameters. 

2) Energy per baryon, charge etc are strongly parameter dependent. 

3) Normally \Z\ <C A, and the ground state strangelets have slightly posi- 
tive charges. 

4) A liquid drop model (Fermi-gas model) explains the general properties. 
Curvature is a decisive effect. 

5) Strangelets have prominent shell structure for low A. 
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FIG. 4. Energy per baryon in MeV as a function of strange quark mass for 
B 1 / 4 = 145MeV. Four sets of curves are shown. For each set of curves the one with 
lowest E I A represents bulk SQM, whereas the others from below are for baryon num- 
bers A = 100, 10, and 6, respectively. Full curves show the ground state strangelets, 
i.e. those with the lowest energy per baryon for the given A. The other curves 
all have fixed charge, Z = 0, but varying strangeness content. Dotted, horizontal 
curves have zero strangeness (S/A = 0, where the strangeness S is minus the num- 
ber of strange quarks). Dashed curves have S/A = —1, and dash-dot curves have 
S/A = —2 (and no rf-quarks). Results like these obtained within the bulk approx- 
imation to the MIT bag model can be scaled to other choices of B by multiplying 
numbers on both axes by B\£. Notice the wide span of possible strangelet masses 
for fixed B, A, Z, and m s , depending on the quark composition of the strangelet! 
The composition of experimentally accessible strangelets is again determined by the 
production mechanism and the decay properties of the excited states; accurate mass 
measurements would be of great importance. 
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6) Many isotopes are almost degenerate in energy. 

7) Decay modes for metastable strangelets are numerous. Lifetimes are not 
predictable at present, but could well be larger (due to Pauli-blocking) 
than those of the strange hadronic matter discussed by Dover elsewhere 
in these Proceedings, thus providing a possible experimental distinction 
(also, strange hadronic matter can not be absolutely stable in contrast 
to strangelets). 

The conclusions above are mainly based on studies within the simplest 
version of the MIT bag model. Effects of zero-point energy and finite a s have 
only been included in a crude fashion, and actual QCD calculations are beyond 
reach. Real strangelets can have non-spherical shape. They are created at a 
high temperature (which tends to increase the energy per baryon). Et cetera. 

Work on some of these complications is in progress, but clearly our theoreti- 
cal understanding of strangelets is still in a rather primitive state. Experimen- 
tal data are desperately needed before constructing even more complicated 
models! 
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